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We study quantum processes with two or more adiabatic parameters. When the param- 
eters are compactified, a derivative of the Hamiltonian is closely related to a first Chern 
number. This topological formulation is applied to the AC Josephson efi'ect and the spin 
Hall effect in semiconductors. 
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1. Introduction 

Recently applications of the Chern number have attracted a lot of interest. The 
notable example is the quantum Hall effect. It can be argued that the precise quan- 
tization of Hall conductance should be related to topological nature of the systems. 
In fact the quantization of the Hall conductance in two-dimensional periodic po- 
tentials was proved by using the Kubo formula by Thouless, Kohmoto, Nightingale, 
and den Nijs^. It was shown later by one of uP^that the quantization is due to the 
topological nature of the problem and Hall conductance is given by Chern num- 
bers in the theory of fiber bundle of differential geometry. The base manifold is T^, 
which is a magnetic Brillouin zone, and the fiber is wavefunctions. Niu, Thouless, 
and Wu proposed an interesting idea of twisted boundary conditions which yield 
Chern numbers even in systems without translation symmetrj^]. Other examples 
include quantum Hall effect in three dimensionJ^^^nomalous quantum Hall effect 
in ferromagnets"^, spin Hall effect in vortex states and Thouless pumping ^. 

In this letter we propose a new method involving Chern numbers to have a 
unified point of view in adiabatic quantum process. We shall apply this theory to 
the examples above. In addition we shall discuss AC Josephson effect and spin Hall 
effect in semiconductor and obtain a new perspective to these problems. 



1 
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2. Formulation 



We start with systems witlr two or more adiabatic parameters. Consider a Hamil- 
tonian H{g{t)) witli A^(> 2) adiabatic parameters g(i). Schrodinger equation is 



d 
i — 
dt 



*g(0)=i^(gW) *g(t))- (1) 



In the adiabatic approximation, the parameters g(t) change in time very slowly. 
Then, the solution of Eq. ([Ij in the adiabatic approximation is 

where |^'ng(t)) is an eigenvector of the snapshot Hamiltonian which satisfies 
Schrodinger equation H{g{t)) \'^ng{t)) — en(g(i)) |^ng(t))- The presence of a gap 
around n-th level and the absence of the degeneracy are assumed. For a later con- 
venience a phase 7n(g(t)) = i /q dt' (^ng(t')l^l^ng(t')) introduced. Note that 
the phases of the states can not be determined uniquely at this stage since |5'ng(t)) 
is the eigenstates of the instantenious snapshot Hamiltonian. When the system 
makes a loop in the parameter space we have H{g{T)) = iJ(g(0)) and eigenfunc- 
tions are same for t = and t = T except for phase difference. This difference is 
called Berry's phasJIH and given by Tnil) = dtjnisit)) = ^fi^S- Ag"\ where 
A.(")(g) = ('I'ngI ^ l^'rag) is a gauge connection in the parameter space '2'. The 
phase 7ra(g(i)) is nonintegrable, it can not be written as a function just of the end 
points, because it depends on the geometry of the path connecting them as well. It 
can be considered as holonomy'^. Namely, Berry's phase is path dependent and 
contains geometrical information as a holonomy. 

Let us consider derivatives of the Hamiltonian with respect to adiabatic pa- 
rameters. Hereafter, we eliminate the unimportant dynamical phases in Eq. ([2|) 
by putting e„(g(t)) — 0. This manipulation is just to redefine the Hamiltonian as 
H{g{t)) ~ e„(g(i)) H{g{t)). This is allowed because eigenstates are unchanged 
and does not affect the argument below. The expectation values of the derivative 
are 



d 

dga 



{(^„g(,)| ii(gw) l^ngit))} - \Hism 



*ng(t)|i?(g(t))| 
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d 

dga 



rig(t) 



dt 



dga 



'Sit) 



9*, 



sit) 



dt 
dga 



ng(t) 



dga 



■sit) 



d^, 



■sit) 



dt 



dt 



sit) 



dga 



dt 



ng(t) 



^ 7 



sit) 



■sit) 



dga 



= -igb{t) 
= -igbit) 



d^nsit) 



dga 



Jig(t) 



jig(t) 



dgb 



dgb 



9* 



ng(i) 



dga 



(3) 



where a, 6 — l,...,Af, repeated indexes are summed over, and the rela- 
tions ((9*„g(t)/(95a I *ng(t)) = - (*ng(t) I 9^'„g(t)/9gQ) and |d«'„g(t)/(it) = 
gb{t) \ d^ nz(t)l dgb) are used. Such a term cannot arise in one parameter systems . 
Note that our formulation is independent of the spatial dimensions or symmetry. 

We consider the cases where the parameters are on the N dimentional tori, 
< <7i < 27r (i = 1, 2, N). Consider N = 2 case. Averaging out the initial values 
of parameters g'*' = g{t = 0) from Eq. ([3]), we see 



/ dHjg) 
\ dga 



where a,b — 1, 2, and 



T2 (27r)2 
C^"'>eabgb{t), 



ng(t) 



dH{git)) 



dga 



■g(*) 



2^2 



(27r)2^ 



Vg X A(")(g) 



27r ■ 



(4) 



(5) 



Thus it is the first Chern number Af^f^^ and quantized in the unit of l/27r. 

The Chern number is renowned as the quantized Hal l condu ctance of the Bloch 
electron in the magnetic field in two dimensions (2D) Here, it should be 

noted that, the derivation of Eq. ([5]) is independent of the real space structure, 
namely, the periodicity and dimensions. The crucial point to derive Eq. ([5]) is the 
presence of the adiabatic parameter g(i) defined on and the validity of taking 
the average on g^*) . Therefore Chern numbers have a potential to appear in a wider 
range in physics. 

The condition for nonzero C*-"-' is as follows. Introduce "parity" transformation; 



H{gi,g2) H{gi,-g2), |^'«(ffi,g2)) |^«(si -S2))> and 



mation; H{gi,g2) H*{-gi,-g2), l^* 



n(si>92 



"(-ffl,-ff2) 



time reversal" transfor- 
These are analo- 



gous to the parity and time reversal transformations in k-space of 2D Bloch electron 
systems, where k is the crystal momentum. To C*-"'' be nonzero, the system should 
not be invariant under these two transformations, since the Chern number counts 
the vorticity of l^'ng) in g-space and the vorticity changes its sign under these two 
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transformations l2lT2l_ xhen, to preserve the nonzero net vorticity, parity and time 
reversal symmetry have to be broken. 

Consider = 3 case. Average out the initial parameters g^'^ = g(i = 0) from 
Eq. ([3]), we obtain 



T3 



(27r)3 



dH{g{t)) 



dga 



"g(t) 



and 



C 



in) 
ab 



d'g 
T3 (27r)3i 



(«)/ 



dga dgb 

where a, 6 = 1, 2, 3. Let us introduce a three- vector, 



(6) 
(7) 



T3 (27r)3ii 

27r 



(2^)^ 



VgxA(")(g)' 
rAAi")(ge), 



where 



2TTi 



VgxA(")(g) 



(8) 



(9) 



is the first Chern number on a two-torus with fixed gc- This is equivalent to Eq. 
([5]). Because of its topological nature, it is independent of gc- Then, we see that 



^ab - l^^^abcMchc ■ 



(10) 



An extension to A?^ > 3 case is straightforward. 

To sum up, in the adiabatic process we have shown definitely the relation be- 
tween the expectation value of the derivative of the Hamiltonian operator and the 
Chern number, namely, Eqs. (j4|) and (O, and Eqs. ([6]) and (fTO]) . 



3. Applications 

We now apply the formulation described above to a number of problems. 



3.1. integer quantum Hall effect in two and three dimensions 

The electric field is given by E = —A. The crystal momentum k moves adiabatically 
as 



L 



(11) 
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Then, the electric current carried by a state with g(t) is 

e dH{g{t)) 

J(gW) - gg ' (12) 

where d denotes the spatial dimensions. For a filled n-th band in 2D, (U) and ([5]) 
give 

Wn = ^J^Ch^vE,., (13) 

i.e., the quantized Hall current in 2e| 1 | 2 | 12 J 
In 3D, we have, from Eqs. (O and ([TU)) . 



= ^^vkE,f^Ch'C\ (14) 

in agreement with Refs. 1^1. These results show that the problem of quantum Hall 
effects in 2D and 3D can be treated in a unified fashion. 



3.2. Anomalous Hall effect 

The relations (jlip and (fT2)l can be used to this system. Therefore, for the filled 
band, we obtain the quantized Hall current from Eqs. (|4]) and ([5]). From Eqs. 
^ and pO]) it is possible to obtain the 3D current ([14]) which is not discussed in 
Ref. 13. 



3.3. twisted boundary conditions 

Niu, Thouless, and Wu pointed out that the Hall conductance of the system with 

twisted boundary conditions (BC) is represented by the Chern number when the 
3 

finite energy gap exists . Here, we show that the same result can be obtained from 
Eqs. ([ll and Compactified parameters < (a,/3) < 27r/L (here, L is the size 
of the system) are introduced to represent the twisted BC Different values of 
parameters shows the different BC. One can see from Ref. ^ that the parameters 
play the role that looks like k of the Bloch electron. Then, one can put 

g(t)/i = (a, /?) - eEt (15) 

(See, Eq. (fTT|) ). The electric current carried by a state with g{t) is also given by 
Eq. (fT2|) . Averaging out g(0) = {aL,PL) is justified since it has been shown that 
the Hall current (conductance) is independent of (a, /?) Then, from Eqs. ([4]) and 
^ the quantized Hall current for the filled n-th level (fT3|) is obtained . The 3D 
extention is also possible. 



3.4. Quantum spin Hall effect in a vortex state 

The quantum spin Hall effect in a vorex state of 2D superconductor and a superfulid 
■^He film has been investigated In these systems the Bogoliubov quasiparticles 
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are in the Bloch state. The driving force for the spin Hall current is not E but Vi?^, 
where is the magnetic field, and we put g{t)/L = k + where k is 

the crystal momentum of the quasiparticle and the Bohr magneton. The factor 
1/2 comes from spin of the quasiparticle. The spin current of the z-component in 
the spin space carried by the quasiparticle is analogous to the charged current of 
Bloch electrons ^ 

, _ 1 dH^dG _ 1 dH^dG 



" 2L2 ak 2L dg ' 

where H^^G jg ^j^g Jjamiltonian of the Bogoliubov-de-Gennes equation in the vortex 
state. Then, from Eqs. ^ and ([5]) we obtain the quantized spin Hall current carried 
by the filled n-th band 

which agrees with Ref. 



3.5. Thouless pumping 

Thouless considered ID Bloch electrons with an adiabatically moving periodic po- 
tential, which looks ac-like external field at fixed position. He showed that dc charge 
transfer occurs and its amount per cycle is represented by the Chern number . The 
same result can be obtained from Eqs. (|4]) and ([5]). We denote L the spatial peri- 
odicity. The crystal momentum k varies < k < 2Tr/L. The Hamiltonian in the 
fc-space has a form Hk{x — Vt), where ^ << 1 is the velocity of moving potential. 
It is obvious that the system is cyclic in time with a period T — L/v. For conve- 
nience, introduce a parameter a varying < a < L and consider a Hamiltonian 
Hk{x + a — Vt) = Hk^a{t) instead of Hk{x — Vt). It is clear that the charge transfer 
during the period T is independent of a. Then, we can average out a. We may put 
g(i) = {Lk,'^{a-Vt)) and Hu^ait) = H{g{t)). By using Eqs. g]) and ©, and 
the current (fT2|) with d = 1, we obtain for a filled n-th band 

(^.)„ = e^. (16) 
The charge transfer per cycle agrees to Ref. ^ 

(AQ)^ ^ r dt{J,)^^eM^l (17) 
Jo 

It is noted that the above formulation seems to have some relation to the discussions 
in Ref. ISI. 

We can clearly see the relatiorP between Thouless pumping and the AC Joseph- 
son effect by using the general formulation we have introduced. Consider two super- 
conductors 1 and 2 coupled by an insulating film. We denote 9i and 6*2 the phase of 
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the order parameter of the superconductors 1 and 2, respectively. The Hamiltonian 
for the relative phase ip = di — 02 is 



Hit) = -4.Ec— - Ej cos(^ - 2eVt), (18) 



where E^ and i?j are the charging energy and couphng energy of the junction, 
respectively fl^. The period of the Hamiltionian T — -k jeV . We suppose that Y << 
1 and use the adiabatic approximation. Since the system is invariant under (p —>■ 
9? + 27r, the eigenstates at fixed t arc in the Bloch states ^g"-* {ip, t) — e^^'^C/^"'' (p, t), 
where < q < 1 and U^''\p,t) = [/^"V + 27r,t). The Haniihonian for u!j"-\(p,t) 
is 

Hg,c,{t) = 4:Ea (^-i^ +<?j - Ej cos{p + a~ 2eVt) 

where a (0 < a < 27r) is an arbitrary parameter. We can average out a since the 
change of the relative phase per cycle is independent of a. From Eq. p9)) . we see the 
one to one correspondence to our formulation. We may put g{t) — {2Trq,a — 2eVt) 
and Hq^a{t) = H{g{t)). We see 



dHisit)) _ 1 5i/,,„(t) _ ip 



(19) 



dgi 271 dq 2tt' 

By using Eqs. ^ and (O, we obtain for the n-th level 

-4^eycf = -2eyA/-(;^ (20) 

The change of relative phase per cycle is 



T 



(A^)„^ / dt{p)^^ = -2^N^-;i. (21) 



This result is formally equivalent with the formula for Thouless pumping Eq. (|17p . 
In our knowledge, the observable physical measurement for Eq. (|2ip is not clarified 
at present. 



4. Spin Hall effect in semiconductors 

It has been aregued that, in semiconductors with spin-orbit couplings, spin current 
may flow perpendicular to the electric field Then this effect is different from 
the spin Hall effect in a vortex state El. Here, we would like to consider a condition 
for a Chern number expression of the spin Hall current in semiconductors by using 
the general formalism obtained in Section 2. 

To discuss the spin Hall effect, we generalize dH{g)/dg 

_UdH{g) \ 
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where Og is an observable and {A, B} — AB + BA. Let us consider the condition to 
have the Chern number expression for the expectation value of generalized velocity 
(vo). When 

[iJ(g),Og]=0, (22) 



we have their simultaneous eigenstates, i.e. H{g) 



^ng^' ), and if the eigenvalue o„g is independent of g, i.e., 



(23) 



we can see that, repeating the derivation of Eqs. ^ and ([5|) for N=2 case with a 
slight modification, 

(-«)n-^E°"^ir"\ (24) 

On 

where A/'^^'^'^'^is the Chcrn number 



AT 



(n,o„) _ f <Pg 



ch 



T2 



27ri 



Vg X A^"'°") 



and A^"^°^^ = ( *ng"^|Vg|*i°g"^ ). Eqs. ^ and ^ are the conditions for the 



Lg — \^ ^ ng I V g I ^ „g ^ 

topological expression of (vq)^. The N(> 3)-parameter extension is possible. 

Consider a semiconductor. We apply the relation (jlip. In case that k is not a 
good quantum number, we apply the twisted boundary condition technique in Ref.l^ 
and use Eq. (fT5|) instead of Eq. pT|) . Put Og = Sz [Sz] ^-component of spin). Then, 
vo is proportional to the spin current ^1^. When the spin-orbit coupling is present, 
Sz does not commute with the Hamiltonian and not conserve, i.e. Eq. (|22p is not 
satisfied. It is possible to define conserved spins with various manipulationJ^QE^, 
but the eigenvalues of them depend on g^, i.e. Eq. is not satisfied. See, 

also Ref. We can obtain the topological expression for the spin Hall current in 
semiconductors, if we succeed to find out a new conserved spin whose eigenvalues 
satisfy Eq. ((23l) . 
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